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This note reports on an independent implementation of calculating one-loop amplitudes 
semi-numerically using generalized unitarity techniques. The algorithm implemented in 
form of a C-|— |- code closely follows the method by Ellis, Giele, Kunszt and Melnikov [T] 
[5]. For the case of gluons, the algorithm is briefly reviewed. Double-precision results 
are presented documenting the accuracy and efficiency of this computation [3]. 

1 Introduction 

An automated next-to-leading order generator for Standard Model processes is highly desir- 
able [4j. With recent developments of generalized unitarity [5] and parametric integration 
methods [5] such a generator seems to be within reach [U [5J [71 [5] • A first crucial step is 
the development of stable and fast algorithms for evaluating one-loop amplitudes through 
generalized unitarity cuts. A C-I-+ code is especially of interest because of the ease with 
which it can be integrated in leading-order generators such as COMIX [3] . The leading-order 
code will then be used to compute the cut graphs. Eventually, such a retrofitted generator 
will be able to generate all necessary amplitudes for a next-to-leading order Monte Carlo 
program for any Standard Model process of interest to the collider experiments. 

2 Construct one-loop by tree-level amplitudes — algorithm in brief 

The full iV-gluon one-loop amplitude can be constructed from the leading colour-ordered 
amplitudes [10], which can be calculated by A^^({pi, Ki}) = A'j^ + Rpf depending on external 
momenta Pi and polarizations Hi. The cut-constructible part reads 

N / J lll2*3U *1*2*3*4 / ^ *1*2*3 *1*2*3 / ^ ni2 ni2 ^ ' 

[«l|«3] [n|''2] 

employing the short-hand notation = 1 < ii < ^2 < • • • < H/ < N , and AI denotes 

the number of cuts. The master integrals are defined as I^'^}.^^^ = J d^£ {iir^/'^ di^ ■ ■ ■ d^jj- 
and the inverse propagators di arc functions of the loop momentum: di{£) = {i + Qi — QimY 
with qk = '^j^iPj- The rational part Rpf is represented by 

^(4) J7) ,2\ 

n \ " "iii2i3i4 , \ ' ^'iii2i3 \ " / [Qil ~ 9*2 j \ 7,(9) /o\ 

- 2^ — Q — + 2^ — — 2^ I Q I &n.2 ' (2) 

[illu] [nl«3] [ilK2] 

cf. [21 [I]- For each possible cut configuration ii ■ ■ -im, the box (^["^21314)' triangle {c'^l^i.^) 
and bubble (^["j^) coefficients appearing above are found as solutions to the parametric form 

* J.W. thanks the Aspen Center of Physics, where parts of this work were accomplished during the Summer 
Workshop "LHC: BSM Signals in a QCD Environment". 
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of the imintegrated ordered one-loop amplitude, 
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dii di2 



(3) 



This decomposition of the integrand A''j!^°' {£) ~ A/''^^' (^) / ((iic?2 • ■ ■ dN) has been generalized 
to higher (integer) dimensionality of the internal particles, i.e. spin-polarization states and 
loop momenta respectively have dimension Ds and D < Ds (as required by dimensional 
regularization). The Ds dependence of the integrand can be eliminated by taking into 
account that the numerator only linearly depends on the spin-space dimension: M^^'^i) = 
Afo{£) + {Ds — A) Afi{i) ■ Moreover, only up to 5-point terms (i.e. M < 5) need to be included 
in the parametrization, since the loop momentum effectively has 4 -f 1 components only: 
A^j^°\£) ~ '4j^°''(£i, . . . , £4, [— i"^]-^/^). These two additions essentially are sufhcicnt 
to disentangle the rational part in the same way as the cut-constructible part [5] . 

The numerator functions of the A/-point terms are polynomials encoding the loop- 
momentum dependence residing in the space orthogonal to the physical space defined by 
the external momenta. The orthogonal space is spanned by the basis vectors n^. The form 
of the polynomials is richer compared to the four-dimensional case. They now include the 
coefficients that determine the rational part, e.g. for the box numerator, one finds 
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where ai 

functions, loop momenta £ ~ £i-^...ij^, 
3 = ii, 



To solve for the coefficients of the numerator 



have to be constructed such that 
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Using the rii vectors, loop momenta fulfilling the constraints can be generated: £ij^...ij^, = 

^ii-iM +J2i=M where = ~^j?--uf -J2i=M+i In general, the solutions are found 
as complex momenta. The vectors 1/^,...^^^ reside in the physical space; they are constructed 
from sums of external momenta as specified by the cuts, cf. [1]. By successively applying 
quintuple, quadruple, triple and double I?s-dimensional unitarity cuts, all coefficients can 
be determined. The cuts yield M onshell propagators factorizing the unintegrated one-loop 
amplitude into M tree-level amplitudes. Hence, the residues can be evaluated: 



Res,; 



where £i^. = 



A^^^'{£) 




(Afc). 



^(Afc-j-i) 



(7) 



lik ^ liM and the sum is over internal polarization states. 



3 CH — h implementation and results 



The algorithm described above has been implemented in a new CH — h code. The only interface 
is to link the QCDLoop package [11] for evaluating the master integrals. The program is ca- 
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pable of calculating the Aj^ {{pi, Ki}) amplitudes in double precision in the four-dimensional 
helicity scheme. This allows for crosschecks with the results obtained in Refs. [7l I12|. 

The construction of the orthonormal sets of the D — M + 1 basis vectors and the Dg — 2 
polarization vectors follows the method outlined in . In addition, the Ui vector generation 
{i = M, . . . ,D) has been set up such that basis vectors obtained for large- Af cuts can be 
re-used for suitable lower-M cuts. The basic strategy, which has been implemented to find 
the coefficients, is as follows: by using the freedom in choosing loop momenta, m solutions 
and, therefore, algebraic equations, such as eq. ([4]), can be generated to solve for n < m 
coefficients^ First the dependence on Dg is eliminated by computing: KeSij^...i^[AN{£)] = 
{D, - 3)ReSi,...,^,[A^]^''\e)] - {D, - i)Res,,...i^,[A]^'^^\e)]. Then higher-point terms are 
subtracted yielding numerator factors ei-i^...i^{£) etc. that are independent of Dg, i.e. eqs. (O 
and (O work without the (Dg) label. For the coefficients of the cut-constructible part, 
one can dispense with the determination of the Dg dependence of the residues and set 
D ~ Dg = 4:. This, in addition, leads to smaller subsystems of equations, which can be solved 
separately, e.g. cq. (jH) simplifies to di^...i^ {£) = d[^\,^^ +ai d[^\,^^ . The tree-level amplitudes 
needed to obtain the residues in cq. ([T]) are calculated with Berends-Giele recursion relations 
[13], adjusted to work for gluons in higher dimensions. For efficiency, currents, which involve 
external gluons only, are stored for re-use in evaluating other residues. 

A number of consistency checks was carried out to verify the correctness of the imple- 
mentation. The gauge invariance of the results and their independence of different choices 
for loop momenta and dimensionalities D and Dg were tested. Coefficients themselves, the 
pole structure of the amplitudes and, finally, the amplitudes themselves have been compared 
to analytic results for various N and different momentum and polarization configurations 
of the gluons. Agreement within the limits of double- precision calculations has been found 
with the numbers produced by Rocket for the fixed phase-space points given in 

In the following, studies are presented that have been conducted to examine the accuracy 
and time dependence of the numerical calculation. 

Accuracy of the results. The quality of the numerical solutions can be estimated by 
analyzing the logarithmic relative deviations, which are defined as 



where the analytically known pole structures of the one-loop amplitudes are taken as ref- 
erence for double (dp) and single (sp) poles, while for finite parts (fp), two independent 
solutions are compared with each other. Figure [T] shows the e distributions together with 
the number of generated phase-space points for various N. The top row of numbers in the 
plots displays the means of the distributions. All results have been obtained for the same 
cuts on external gluons as reported in 0; for the effect of tighter cuts, see Figure [3] (right). 

For N = 15, the double-precision evaluation of the coefficients clearly is not sufficient to 
yield reliable finite-part result. The loss of precision as N increases is correlated with the 
more frequent appearance of small denominators and large numbers characteristic for the 
calculation. The two rightmost bottom panels of Figure [T] present two examples by depicting 
the range of magnitude taken by Gram determinants, used to evaluate the T^j...^^,^ vectors 
of external gluons, and e^^.l.jg coefficients, which can be of 0(10^^) for N = 15. Double 

''In principle, an infinite number of equations can bo generated to fit a fixed number of unknowns. 
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Figure 1: Double-, single-pole and finite-part accuracy (in double precision) of the +H . . . 

one-loop amplitudes for = 6, . . . , 11, 15 gluons; see also text and right panel of Figure [3l 
Bottom row: center, double-logarithmic distributions of Gram determinants involving sets 
of external gluons and, right, ef\,^ coefficients for the TV = 6, 8, 15 gluon setups. 



precision will then be insufficient to make cancellations as they may occur e.g. in eq. ([6]) 
manifest. The scatter plots of Figure [2] visualize that (partial) correlations exist between 
the relative accuracy of the finite part and the smallest Gram determinant of extcrnal- 
gluon sets, the largest e\^\.i^ coefficient and the single-pole accuracy. In all cases, the areas 
of scatters shift with increasing N towards worse accuracy and more extreme values of 
Gram determinants and 5-point coefficients. The calculation may still involve other small 
denominators, such as the leftover dj in the subtraction terms of e.g. eq. This leads to 
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Figure 2: Finite-part accuracy versus minimal Gram determinant of external-gluon sets 
(left), maximal e|^'..j^ coefficient (center) and single-pole accuracy (all in double precision) 
for the = 6, 8, 11 gluon setups of above. 

instabilities (even for small coefficients) and uncorrelated areas in the plots are populated. 
Given these correlations, it can be seen that one way of achieving higher accuracy is to 
compute the coefficients in quadrupole precision. This has been pointed out in Ref. 0. 

Efficiency of the calculation. As estimated in Ref. [3] (p. 31), the algorithm is expected 
to have polynomial complexity, see also [7] . The computing time r^r to calculate an ordered 
A^-gluon one- loop amplitude should scale as A'^^ -|- . . . , with the leading term having x = 9 
dominating the behaviour for large N. The results for rjv with = 4, . . . , 20 arc shown in 
Figure [3] together with the exponents xn ~ In ^ / In . The plots demonstrate that 
the CH — h algorithm as implemented displays the predicted scaling. 
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Figure 3: N dependence of the computing time; iTv-cxponents (center), see text. Times 
refer to using a 2.20 GHz Intel Corc2 Duo processor. Tighter gluon cuts were used: \rii\ < 2, 

iij 



P±.i s > 0.1, Rij > 0.7, denoted as in 0. The last plot shows the accuracy improvement. 



4 Conclusions 

It has been shown that the generalized unitarity method of Refs. [HIl] can be implemented in 
a stable and fast C++ program. The one-loop A^-gluon amplitude has been used as a testing 
ground. The next step is to integrate the C++ code in an existing leading-order generator. 
The resulting upgraded generator will be able to generate both virtual and bremsstrahlung 
contributions for arbitrarily complex Standard Model processes. The final step towards a 
full next-to-leading order Monte Carlo is to add the necessary phase-space integrations. As 
was shown in Ref. |14j the virtual matrix elements calculated with the generalized unitarity 
method can be used in ncxt-to-lcading order Monte Carlo programs. 
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